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Abstract 

We present improved upper bounds for the size of relative {p, £)-approximation for range spaces with 
the following property: For any (finite) range space projected onto (that is, restricted to) a ground set of 
size 71 and for any parameter 1 < k < n, the number of ranges of size at most k is only nearly-linear 
in n and polynomial in k. Such range spaces are called "well behaved". Our bound is an improvement 

over the bound O ^ ^ inti'oduced by Li et al. [17] for the general case (where this bound has 

been shown to be tight in the worst case), when p e. We also show that such small size relative 
(p, e)-approximations can be constructed in expected polynomial time. 

Our bound also has an interesting interpretation in the context of "p-nets": As observed by Har- 
Peled and Sharir [13], p-nets are special cases of relative {p, £)-approximations. Specifically, when e is a 
constant smaller than 1, the analysis in [13, 17] implies that there are p-nets of size 0(log {l/p)/p) that 
are also relative approximations. In this context our construction significantly improves this bound for 
well-behaved range spaces. Despite the progress in the theory of p-nets and the existence of improved 
bounds corresponding to the cases that we study, these bounds do not necessarily guarantee a bounded 
relative error 

Lastly, we present several geometric scenarios of well-behaved range spaces, and show the resulting 
bound for each of these cases obtained as a consequence of our analysis. In particular, when e is a 
constant smaller than 1, our bound for points and axis-parallel boxes in two and three dimensions, as 
well as points and "fat" triangles in the plane, matches the optimal bound forp-nets introduced in [3, 25]. 
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1 Introduction 



Motivated by the problem of approximate range counting, relative (p, e)-approximations have been intro- 
duced by Har-Peled and Sharir [13], where they revisited the seminal work of Li et al. [17], and showed how 
to apply it in order to derive known bounds on the size of various notions of sampling, and, in particular, 
relative (p, e) -approximations. We recall their definition: A range space (It, 3?) is a pair consisting of an 
underlying universe U of objects and a certain collection 3? of subsets {ranges). Of particular interest are 
range spaces of finite VC-dimension; the reader is referred to [15] for the exact definition. Informally, it 
suffices to require that, for any finite subset X C U, the number of distinct sets r n 31 for r £ "Rhe 0{\X\'^) 
for some constant d (which is upper-bounded by the VC-dimension). This is indeed the case in many geo- 
metric applications. In a typical geometric setting X is a subset of some infinite ground set U, e.g., U = M'^, 
for some low or fixed dimension d, and thus X is a finite point set in W^, and 31 is a collection of subsets 
(ranges) of U of some simple shape, such as halfspaces, simplices, balls, ellipsoids, and boxes in W^, where 
d is considered a fixed. In general, range spaces involving semi-algebraic ranges of constant description 
complexity, i.e., semi-algebraic sets defined as a Boolean combination of a constant number of polynomial 
equations and inequalities of constant maximum degree, have finite VC-dimension (see, e.g., [12] for further 
details and examples). In what follows we assume to have a finite range space defined on a set X of objects, 
and, with a slight abuse of notation, denote its set of ranges by 31. 

Following the notation in [13], the measure of a range r G 31 is the quantity' = Given two 

parameters, < p < 1 and < e < 1, we say that a subset Z C X is a relative {p, e) -approximation if it 
satisfies, for each range r, 

X(T)(l-e) <Z{t) <X{T){l + e), if X(t) > p, and 

X{t) — ep < Z{t) < X{t) + ep, otherwise. 

In fact, a slightly more general notion is the so-called {v, a)-sample [12, 14, 17], in which case the subset 
Z Q X satisfies, for each range r, 

jZ(r)-X(r)| 
Z{t) + X{t) + z/ 

As observed by Har-Peled and Sharir [13], relative (p, e)-approximations and {v, a)-samples are equivalent 
with an appropriate relation between p, e, and v, a (roughly speaking, they are equivalent up to some 
constant factor). Due to this observation they conclude that the analysis of Li et al. [17] (that shows a bound 
on the size of (z/, a)-samples) implies that, for range spaces of finite VC-dimension d, there exist relative 
(p, e) -approximations of size ^^^^^p^p^, where c > is an absolute constant. In fact, any random sample 
of these many elements of X is a relative (p, e)-approximation with constant probability. More specifically, 
success with probability at least 1 — g is guaranteed if one samples '^^'°g(i/p)+^°g(i/'?) elements of X. These 
bounds are in fact an improvement over the bound '^'^^"g i^/(p^))+^°s (i/q) obtained by Haussler [14] and 
Pollard [26]. ^ ^ 

It was also observed in [13] that p-nets and e-approximations are special cases of {u, a)-samples. The 
first is a subset N X with the property that any range r G 3i with |r n X| > p|X| contains an element of 
N, and the latter is a subset Z C X with the property that any range r G 31 satisfies: \X{t) — Z{t)\ < e. 



'Technically, t C X , so \X f] t\ = |r|. However, (i) we also use this notations for subsets Z C X, and (ii) typically X is a 
finite subset of R'* and the ranges are best described as intersections of simply-shaped regions with X, as above. 
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In this paper we present improved bounds on the size of relative (p, e) -approximations under certain 
assumptions, and emphasize their implications to p-nets — see below. 

Previous results. It has been shown in [29] that range spaces of finite VC dimension d always admit an 
absolute-error e-approximation of size 0(log {l/e)/e^), where the constant of proportionality depends on d 
(see also [5, 6, 20, 24]). In fact, a random sample of that size is an e-approximation with constant probability. 
As noted in [13], the analysis in [17] (see also [28]) reduces this bound to 0(l/e^), where, once again, a 
random sample of that size is an e-approximation with a constant probability. In fact, as shown in [5, 6, 20], 
the size of the e-approximation can be further improved to be slightly higher than 0(l/e^~^/('^+^)); see [13] 
for a more comprehensive review of these results. 

Concerning relative {p, e)-approximations, as noted in [13], any absolute error ep-approximation Z will 
approximate ranges of measure at least p to within relative eixor e. Nevertheless, the bound in [29] for 
absolute-enw approximations just mentioned yields a sample of size 0(log (l/(ep))/e^p^) in this case, 
whereas the bound of Li et al. [17] is smaller by roughly a factor of 1 /p (see the discussion above). 

In addition to the observations made by Har-Peled and Sharir [13], their analysis also improves the 
bound on the size of relative {p, e)-approximations for several special cases of geometric range spaces. 
Specifically, the bound obtained for point sets in the plane and halfplane ranges has been improved to 
Q ^ log /^^(i/(p£)) ^ improvement if e is sufficiently small with respect to p). In 3-space they 

reduced the bound to O ^ ^"g ^^Wip^^)) ^ ^ although this latter case is somewhat restricted, as this is the bound 

on the overall size of 0(log subsets, and each half space range has one such subset that constitutes its 

relative approximation. 

Our results. In this paper we derive improved bounds for range spaces with certain favorable properties, 
which we call "well-behaved". Our goal is to improve the dependency on p, which is somewhat different 
than the improvement obtained in [13]. That is, for any induced (finite) ranges space, the number of ranges 
of size k is only nearly-linear in the size of the space and polynomial in k, where > is an integer 
parameter. We first present our technique on abstract range spaces that satisfy these properties, in which 
case we show there exists a sample of size roughly 0((loglog + log (l/e))/e^p) (see Corollary 2.5 
for the exact bound), from which the relative approximation is derived. We also present an algorithm that 
constructs such a sample in expected polynomial running time. In fact, this sample consists of two subsets 
of the aforementioned overall size, on which we assign weights, in order to obtain a single weighted sample 
with the desired properties; see Section 2 for this derivation and the discussion below. 

As observed by Har-Peled [12] relative approximations are interesting in the case where p = o{e), since 
one can approximate ranges of measure larger than p with a sample that has only linear dependency on 1/p. 
Our bounds improve over the bound of Li et al. [17] for these cases. Nevertheless, our bounds also have 
an interesting interpretation when e is some fixed constant in (0, 1). In this case, as mentioned earlier, the 
relative (p, e)-approximation becomes a p-net. Thus when e is a constant smaller than 1, a p-net of size 
0(log (l/p)/p) (with a constant of proportionality depending on the VC dimension d) always exists by 
applying the bound of Li et al. [17]. In this case, each range of measure at least p contains roughly log (1/p) 
points of the sample. Nevertheless, in some cases it might be wasteful to sample this number of points in 
order to guarantee a small relative error. There are several known constructions for small-size p-nets, see, 
e.g., [3, 8, 19, 16, 21, 27, 30] for several such improved bounds. However, we are not aware of any such 
small-size p-nets that are also relative approximations. Our bound guarantees such properties. For example. 
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for primal range spaces of points and axis-parallel boxes in two and three dimensions this bound becomes 
O ^°g(^/p) ^ instead of the standard bound 0(log {l/p)/p). This bound matches the previously known 
(optimal) bound for p-nets [3, 25]. 

We also note that the task of constructing a relative (p, e) -approximation is much more difficult than 
just constructing p-nets. Indeed, in p-nets we aim to add sufficiently many objects (e.g., points) into the 
output sample so as to guarantee that each range of measure at least p is indeed stabbed, whereas in relative 
(p, e) -approximations we also need to keep the set of the chosen objects under control in order to have a 
balance between the original measure and the approximate measure for each range r. 

In our analysis we initially replace the set of the input objects with a relative (p, e)-approximation 3" of 
size O ■ We then classify the objects of 3" as "heavy" (that is, objects that participate in many 

ranges) and "light" (otherwise). We show that the number of heavy objects is relatively small. For the light 
objects, we randomly and independently choose each of them into the new relative (p, e) -approximation 
with some probability vr (calibrated to produce a sample of a certain expected size). Since the original range 
space is well-behaved, we are able to show that each range r (of any measure < < 1) admits a small 
degree of dependency in the space of events ''Ar'- does the relative error for r exceeds eT. We then apply 
the asymmetric version of the Local Lemma of Lovasz in order to conclude that with a positive probability 
there exists such a "good" sample for which the relative error of all such ranges does not exceed e. In fact, 
we show that both conditions of having a good sample (in the above sense) and keeping its size close to 
its expectation (up to some constant factor) happen with a positive probability — this latter condition can be 
included into the Local Lemma, which extends to that case. See Section 2 for these details. 

The manner in which we construct our relative approximation enforces a somewhat different form for 
the measure Z{t). This is due to two main ingredients of our analysis: (i) Contrary to the standard construc- 
tions of relative {p, e)-approximations [17] (as well the the construction in [13]), absolute e-approximations 
(e.g., [5]), and p-nets (e.g., [8, 30]), where the size of the sample is fixed^, in our probabilistic model each 
object is chosen independently with a fixed probability vr, and thus in the denominator of the measure we 
replace the size of the sample by its expectation, (ii) Our sample consists of two subsamples, where the first 
one H consists of a pre-determined subset of the input, which, due to its small size, is taken in its entirety 
into the output, and the second sample 3"i is obtained by choosing each input object randomly and indepen- 
dently with probability vr . This results in a biased sample. Specifically, we obtain a weighted measure, in 
which we assign a unit weight to each object in 3"i , and a weight vr to each object in H. 

Lastly, we list several useful applications of well-behaved geometric range spaces, including primal 
range spaces of points and halfspace ranges in two and three dimensions, points and axis-parallel boxes in 
two and three dimensions (as well as points and "fat" triangles in the plane), and dual range spaces involving 
planar regions of nearly-linear union complexity. 

We note that our technique is inspired by the machinery of Varadarajan [30] for constructing small-size 
p-nets for dual range spaces of "a-fat" triangles and points in the plane, where the idea for exploiting the 
(simpler version of the) Local Lemma of Lovasz has initially been introduced, as well as classifying each 
object as "heavy" or "light". Nevertheless, the technique in [30] does not necessarily produce a relative 
approximation, but only guarantees that the sample is a p-net, which is the reason we had to generalize and 
enhance the ideas in [30] in order to be matched with the more intricate scenario arising in our problem. 

^In the sampling model in [17] the size of the sample is always bounded by O ^^r^^j but the objects are chosen with 
repetitions. 
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2 The Construction for Weil-Behaved Range Spaces 



2.1 Preliminaries 

Well-behaved range-spaces. Let (It, 3?ix) be a range space of finite VC dimension. We say it is well 
behaved it it has the following property: Let {X, 3?) be any range space projected onto a finite subset 
X C U, where 3? = {{r n X} \ r G 3?u})' ^rid Put n := \X\. Then, for any pai^ameter 1 < k < n, the 
number of ranges in 01 of size < A; is at most 0{n(f){n)k^), where 0(-) is a slowly-growing function, and 
c > is an absolute constant. In other words, for any induced (finite) range space (X, 3?) the number of 
ranges of size at most k is only nearly-linear in |X| and polynomial in k. Note that by definition any induced 
finite range space {X, 3?) as above is also well -behaved^ . 

In what follows we assume, without loss of generality, that < p < 1/8. Otherwise, if we also have 
e > 1/8, then the size of the relative approximation is a constant, and if0<e<l/8is arbitrary (and 
p > 1/8) then an (absolute) (e/8)-approximation always yields an error smaller than ep. In this case the 
size of the sample is only 0(l/e^), as shown by Li et al. [17]. In addition, we assume p < e, otheiAvise, we 
output a sample of size 0(log {\/p) / e^p), as shown in [17]. 

2.2 The Construction 

Let {X, 01) be a well-behaved (finite) range space. We first replace the objects in X by a sample 3" C X that 
is a relative {p, e) -approximation for (X, 3?); let The the resulting collection of ranges projected onto 3". By 
the discussion in Section 1, there exists such a sample of size '^^"^ia^'^^''^ where I? > is a sufficiently 
large absolute constant. Moreover, a random sample of X of that size is a relative {p, e) -approximation with 
constant probability.^ The replacement of X with 3" implies that each range r G 3J satisfies: 

' ^'l-e)<^——^<^——^{l + e) ifX{T)>p, 



\x\ ^ ' - I j| - \x 

and 



< — fTPl — < — j — j h otherwise. 



\x\ " - I j| - |X| 

From now on, we focus on the construction of an improved relative (p, e)-approximation for (3^, T). 
This (standard) reduction involves no loss of generality, because, as is easily verified, a relative (p(l — e), e)- 
approximation of a relative (p, e) -approximation is a relative (p, 2e + e^) -approximation. Hence by scaling 
p and e by the appropriate constant factors, it suffices to construct a relative (p, e) -approximation for (3", T). 

Our construction partitions r n 3" into two subsets, and represents the approximation for |r n 'J\ as 
the sum of two appropriate sample measures — see below for details, and for the consolidation of the two 
samples into a common (weighted) set. 

A range r is said to lie at the ith layer Sj of T, for z = 1, . . . , log (1/p), if it satisfies 

2-ip < < 2>. (1) 



^We note that just the fact that (11, %\x) is well-behaved already implies that it has a finite VC dimension. Indeed, by definition 
we have that the total number of ranges in any projected range space (X, K) is polynomial in This implies that the so-called 
shattering dimension is finite, and thus the VC dimension is also finite. See [12] for these straightforward details. 

■*Note that this bound is somewhat suboptimal, as the log (1/e) factor in the enumerator can be removed by the analysis of [17]. 
However, due to technical reasons imposed by our analysis, we need a slightly larger sample, which is the size that we choose. 
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For the sake of completeness, the 0th layer Sq consists of those ranges r with ^-^j^ < P- 

Consider a fixed layer Sj, where i > 1. Put Aj := — — l^^iiZiEii) Equation (1) and the bound on |F| 
then imply that each range r G Sj satisfies Aj < |r n 3"| < 2Aj. If r G So> then Equation (1) implies 

< |rn3"| < Ai = 2Ao. 

In other words, all ranges r G Sj have size at most 2Aj (and at least Aj, if i > 1), i = 0, . . . , log (l/p), 
and they appear now as ranges from the range space {X, 31) projected onto (that is, restricted to) 3". The 
assumption that the (original) range space {X, 01) is well-behaved implies that the number of the distinct 
ranges in a fixed layer Sj is only OdJjc/^dJDAj'^). 

A classification of the objects. The discussion above implies that the number of containments between 
the objects in 3^ and the ranges in layer Sj is 0(|3"|(^(|9"|)A,^+i), for i = 0, . . . , log (l/p). We say that an 
object in 3" is heavy in Sj if it appears in at least A ■ (/)(|3"|)Aj'^"'"^ ranges (at layer Sj), for a sufficiently large 
constant A> that depends on the constant of proportionality in the bound on the number of containments. 
Otherwise, this object is said to be light in Sj. An easy variant of Markov's inequality implies that the 
number of heavy objects (in Si) is at most 0(|9"|/Aj) = O ^ 2'-^ip ) • Hence, the overall number of heavy 
objects, over all layers i = 0, . . . , log (l/p), is only 0{l/p). Let H denote this subset, which, from now on, 
we just refer to as the heavy objects. Note that, by construction, each remaining object (an object in 3" \ H) 
is light in each layer i = 0, . . . , log (l/p). Put L := 3^\H and refer to its elements as the light objects. 

We next consider, for each range r, the two subsets n r, L n r of heavy and light objects, respectively, 
which r stabs, and approximate each of their measures in turn. First, since the number of heavy objects is 
only 0(l/p), we put all of them in the target relative (p, e)-approximation. In the sequel we describe our 
approximation with respect to the light objects. 

2.2.1 The Analysis for the Light Objects 

We first observe that the analysis above regarding the size of \H\ allows us to assume that \L\ > |3"|/2. 
Indeed, this easily follows by choosing the constant A sufficiently large. 

We now restrict the range space (3", T) to L. Our goal is now to approximate L{t). We keep associating 
these ranges with the same layers Sj, as defined in (1). Note that a range r at layer i may now satisfy 
|t n L| <^ Aj. This can happen, for example, when the overall majority of objects in r n 3" are heavy 
and then |r n L| is considerably smaller than |r n 3"|. These differences, however, do not affect the final 
approximation for F{t), which is a key observation in our analysis; see Section 2.3 for this derivation. 

We next sample each object in L independently with probability vr := ™^^{i°g ^^^^^ , 

and let be the resulting sample. Thus its expected size is O ( max{iogiog(i/p),iog0(|F|)}+iog(i/e) y ^^^^ 



that, by assumption 0(-) is a sublinear function, and p < e, thus we always have vr < 1, as is easily verified. 
The main ingredients of the analysis are shown in the following proposition: 

Proposition 2.1. With some positive probability, 3"i satisfies, for every layer Sj and for every range r o/Si, 
z = 1, . . . , log (l/p), (where Exp(-) denotes expectation)^ : 



^Note that due to our sampling model we replace |?i | by ExpjJ'i}, which is in fact tt ■ \L\. 



5 



and 



knL| ^ ^ . 1 ^ ^ irnJil ^ |rnL| _ ^ . ^ 



e • (2'~ p) < ttpV ^ iri + e • (2*^ P), otherwise. 



\L\ ^ - Exp{Ji} - |L| 

When i = Owe have for each r G So (in which case |r n L| < p|3"| j; 

IrnLl IrnJil IrnLl 



\L\ - Exp{3-i} - \L\ 

Proof. Fix a layer Sj and a range r € Sj, i > 1. Let A^- be the event that 3"i does not satisfy (2) for r. We 
consider separately the following two cases: |t n L| > 2*~^p|3"|, and |t n L| < 2*~^p|3~|. 

(i) |r n L| > 2*-V|J|: Since Exp{|Ji|} = tt\L\, the event Ar in this case IS 

|rn3"i| , |Tn3"i| -, , 

' '<L(T)(l-e) or J — — ^ > L(r)(l + e). 



vrlLl ^ ^ 7r|L| 

That is, IrnJil < 7r|rnL|(l -e) or ImJi] > 7r|rnL|(l + e). Since Exp{|rn Ji|} = vr • |rnL|, 
we have: 

Prob{yi^} = Prob{|rn Jil < Exp{|r n Ji|}(l - e)} + Prob{|r n Ji| > Exp{|r n S^ilKl + e)}. 

Using Chemoff 's bound (see, e.g., [1]) and the fact that |t n L| > |3"| • 2*^^p, we thus obtain: 

Prob{yi^} < 2cxp{-£2Exp{|Tn Ji|}/3} = 2 cxp {-tteV n i|/3} < 2cxp{-7re2 • 2'"V|J|/3} . 

Substituting I J| = 2}2E^M)1, ^ = max|loglog(l/p),l^^^^^^ ^^^^^.^ 

Prob{yi^}<2exp<'-^-2*-i-(max{loglog(l/p),log(/)(|J|)} + log(l/e))[> < ' 



log(l/p)0(|J|) 



where i? > is a constant that depends linearly on D, and can be made arbitrarily lai^ge by choosing D 
sufficiently large. 

(ii) |r n L| < T-~^p\J\: Here Ar is the event: 

^-r,, _ knJil , , .1 



< L(r) - 2*"^ejj or ^ — — ^ > L(r) + 2*"^ep. 



ttILI ^ ^ ^ ttILI 

That is, |rn Jil < 7r|TnL| - 2'-^ep7r|L| or |rn3"i| > 7r|rnL| + 2'-^ep7r|L|. We then have: 
Prob{yi^} = Prob{|rn Jil < Exp{|r n Ji|} - 2'-i£p7r|L|} +Prob{|Tn Ji | > Exp{|Tn Ji |} + 2'-iep7r|L| } 

= Prob||rn Jil <Exp{|Tn (^l-^^-^j^^|+Prob||rn Jil >Exp{|rn Ji|} (^1 ' '^^'^'^ 



|rnL| ; 



Applying once again Chemoff 's bound, we obtain: 
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Prob{.l.} < 2exp - (^-^ j • vr|. n L|/3 = 2exp | 3^ n L| | " 



Since in tiiis case |L n t| < 2*-V • |3"|, and, as stated above, \L\ > |3"|/2, we obtain that the latter term is 
bounded by 2 exp |— vr • ■ (2*~^p) • |3~|/12}, on which we can derive the same bound as in case (i), using 
similar considerations. Hence in summary we obtain in both cases 

The case i = (or rather r G So) follows by similar considerations as in case (ii) above when replacing 
2*~ip with p. This yields the bound Probjyir} < (^ iog(i/p)0(|F | ) ) ^^^^ 

Applying the Asymmetric Local Lemma of Lovasz. We next apply the Local Lemma of Lovasz (see, 
e.g., [1]), to show that 

log (1/p) 

/\ /\ (1 - Prob{yi,}) > 0. 

i=0 T6Si 

This will imply that there exists a sample 3"i of L that approximates |t n I/| as in (2), for all ranges r. 
Specifically, we are going to apply the asymmetric version of the Local Lemma, stated below in the context 
of our problem. 

We first observe that for a pair of ranges r, r', the corresponding events Ar, Ar' are mutually indepen- 
dent if and only if there is no object in L that participates in both r, t'. Indeed, since we sample each object 
of L independently, the two corresponding events Ar, Ar' can affect each other only if there is an object in 
L that r, r' share. In what follows we denote a pair E, E' of mutually dependent events by iiJ ~ E'. 

Let £i denote the collection of events Ar for all ranges r at a fixed layer i, and let £ denote the entire 
collection IJ!^o^^^^^ order to apply the asymmetric version of the Local Lemma we need to show there 
exists an assignment x : £ — )• (0, 1), such that 

Prob{yi,} < x(yi,) • J] (l-x(^,0), (4) 

I'^-j-f 

for each Ar G £. The Local Lemma of Lovasz then implies that 

/\ (1 - prob{>i,}) > n (1 - ^('^-)) > 0- 

Ares. ArGS. 

(Once again, see [1] for further details.) In Lemma 2.2 we show that there exists such a valid assignment. 
This will complete the proof of the proposition. □ 

Lemma 2.2. The assignment :x.{Ar) = exp {2'+^} • Problyi,-}, for each t ^ §i and for each layer Sj, 
« = 0, . . . , log (1/p), satisfies (4) for every Ar G £. 

We postpone the proof of Lemma 2.2 to Appendix A.l, but leave several remarks below. 
Remarks: 1). We note that for each r G §j, the exponent 2^^^ in the bound on Probj^lT-} "beats" the term 



"^"""^ in the degree of dependency (see Appendix A.l). Nevertheless, when i > j, this exponent cannot beat 



Aj, which is the reason we set x{Ar) = exp 2*+ ProblyL,-}. This also demonstrates the crucial property 
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of classifying the ranges according to their layers, and then bounding the probability to fail to produce a 
relative {p, e)-approximation in each of these layers, as we did in Proposition 2.1. In other words, just the 
information |r n L| > or |r n L| < (which is the standard "cut-off" in relative approximations) is 
insufficient to produce a relative bounded eiTor in the manner that we do. 

2). We note that applying the simpler version of the Local Lemma of Lovasz, in each fixed layer i, is almost 
immediate (a similar step has been taken in [30] for "fat" triangles in the plane and points). Indeed, each 
range has a degree-of-dependency 6 that is at most 0((/)(|3"|)A/^~'^). Following Inequality (3), we obtain 
6 • ProbjyiT-} < 1/e, for a sufficiently lai^ge choice of D (and thus of B). In this case, the Local Lemma 
implies that, with a positive probability, aU the complementary events Ar (for Ar in layer i) are satisfied. 
Nevertheless, this property is not guaranteed for the entire set of events over all layers, as the interaction 
among events from different layers may involve a higher degree of dependency. This is the main reason we 
had to resort to the asymmetric version of the Local Lemma, which, as our analysis shows, overcomes this 
difficulty, and eventually obtains a single sample for all layers. 

Bounding the size of the sample. As noted above, the expected size of 3"i is O ^ ™^^{^°g ^^^^^'ji^ 0(1 J"l}+ipg (lA) ^ _ 
Nevertheless, we need to show that a sample of that actual size exists and that it satisfies the assertions in 
Proposition 2.1. In Appendix A.2 we show that Lemma 2.2 can be extended to include the event 23 that |3"i| 
deviates from its expectation by some constant factor. This yields: 

Corollary 2.3. There exists a sample 3'i L that satisfies the assertion of Proposition 2. 1, whose size is 

Q ( maxjlog log (1/p) ,tog )}+log (1/e) \ 



2.3 Deriving the Relative Approximation 

We now combine the two samples and H that we have constructed in order to derive the relative approx- 
imation for each range r. By construction, |Tn3"| = |TnL| + |Tni7|. Combining the two cases in (2) for 
each range r in Sj, z = 1, . . . , log and adding the term |r fl for each side of the inequality, we 

have: 

|rnL| ^ ^ \Tr\H\ |rn3"i| \Tr\H\ \rr\L\ i^i-i ^ . \Tr\H\ 



■ .{l-e)-e-{2'-'p)+'——^ < ^ ^ + L^^ < .(^l + e) + e■{2'-'p) + ■ 



\L\ ' ' ' \L\ - Exp{Ji} |L| - |L| ^ ' ^ ' ^ ^1 ' 1^1 ■ 

Substituting Exp{9"i} = vr|L|, we obtain: 

|r n 3"! - |r n L|e - e • (2^-ip)|L| < V ^'^^^ r\R\ < |^ ^ ^| ^ |^ ^ ^ ^ . (2^-ip)|L|. 

vr 

Since L C 3", the above inequality can be written as 

Irn J| -elrnJl (2^-^)1^! < 1^ ^ + ^ -^1 < |r n J| + ejr n + e • (2^-ip)| J|, 

vr 

or 

- e) - e . (2-V) < < ^(1 + e) + e . (2-^). (6) 

For i = 1, . . . , log (1/p), since ^pj^ > (2*~^p) (by definition), this implies that 
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When i = 0, we have ^^j^p < ^"^^ then, using similar considerations as above, the case i = of (2) 
implies: 

|rn3"| IrnJil + 7r|Tni?| IrnJI 

— ^\ — 

and we can replace the tenii ep by 2ep in order to be consistent with the previous form obtained for i = 
l,...,log(l/p). 

Note that the measure approximating ^^j^p can be interpreted to be defined on a weighted sample, 
where each object of !Fi is assigned a unit weight, and each object of H is assigned a (fractional) weight 
TT. Also, observe that the total expected weight of the sample satisfies Exp{|3~i| +7r|//|} = 7r|L| + 
7r\H\ = 7r\3'\ (matching the denominator in our measure), and so in our construction this weighted measure 
replaces the standard "uniform measure" Z{t) (defined in the introduction) resulting when the entire relative 
approximation is obtained as a unifomi sample. 

Scaling the parameter e appropriately, we conclude: 

Theorem 2.4. Let (3", T) be a well-behaved range space. Then there exist two subsets 9"i C 3", C 3", 
with the following properties: Each range r G 7 satisfies: 

\t n J| _ ^ Irn Jil + ^lrnffl ^ |rn ^ 



if ^ P' '^'^'^ 



\3-\ ' '- vrlJl - |J| 



|rn3"| Irn Jil +7r|Tni?| IrnJI 

ep < — < — — h ep, 



vrlJl - |J| 

otherwise. Thus the sample !JiU H is a (weighted) relative (p, e) -approximation for (3", T), and its overall 
size is only 0((max{log log log (/)(|3"|)} + log {l/e))/e^p). 



Scaling appropriately the parameter e once again and also p (see Section 2 for the discussion) we obtain: 

Corollary 2.5. Let (X, Jl) be a well-behaved range space. Then there exists a (weighted) relative a (p, e)- 
approximation for [X, JV), of size 0((max{log log (1/p), log (?l)(|3~|)} + log (l/e))/e^p). 



A polynomial-time algorithm. We note that classifying the objects as light or heavy can be done in time 
polynomial in |3"| (we omit these sti'aightforward details here). In order to apply the Local Lemma in a 
constructive manner, we resort to a recent result by Moser and Tardos [23]. In our scenario we apply the ex- 
tended version of Lemma 2.2 (including the event 23), discussed briefly above and proved in Appendix A.2. 
The main property required in order to apply the randomized algorithm described in [23] is the fact that the 
objects in 3"i are chosen randomly and independently with probability vr. That is, the set L induces a finite 
set of mutually and independent random variables. Then each event Ar € £, as well as S, is determined by 
these variables. Omitting any further details, we obtain: 

Theorem 2.6. Given a well-behaved range space (X, one can construct in expected polynomial time, 

a (weighted) relative {p,e) -approximation for {X,Jl), whose size is 0((max{log log (1/p), log (/)(|3"|)} + 
log{l/e))/e^p). 



9 



Applications Using the bound in Theorem 2.4 we obtain several geometric settings that admit small-size 
relative (p, e) approximations. We review these settings and their analysis in Appendix A.3, and conclude: 

Corollary 2.7. Any range space of points and axis-parallel rectangles in the plane admits a ( weighted) 
relative {p, e) -approximation of size 0((log log (1/p) + log {I /e))/e'^p), for any 0<p<l, 0<e<l, 
which can be constructed in expected polynomial time. The same asymptotic bound holds for points and 
axis-parallel boxes in three dimensions, and points and a-fat triangles in the plane (where a > is a 
constant). When e is a constant these bounds match the optimal ©(log log (1/p) /p) bound for p-nets. 

Corollary 2.8. Any range space of points and half spaces in two and three dimensions admits a (weighted) 
relative [p, e) -approximation of size 0((log log (1/p) + log {l/e))/e'^p), for any 0<p<l, 0<e<l, 
which can be constructed in expected polynomial time. 

Corollary 2.9. Any dual range space defined on ( i) pseudo-disks, ( ii) a-fat triangles, ( Hi) locally 'y-fat ob- 
jects and points in the plane admits a relative (p, e) -approximation of size 0((log log (l/p)+log (l/e))/e^p), 
for any 0<p<l, 0<e<l, which can be constructed in expected polynomial time. 

3 Concluding Remarks. 

This study raises several open problems and further improvements, some of which are under on-going 
research. First, it is very likely that the sampling scheme that we introduce can be applied over iter- 
ations, where at the A;th iterations we are given a sample 3^k~i, from which we extract the sets H^-i, 
Lfc„i of the con^esponding heavy and light objects, and then sample each object from Lk_i with probabil- 

ity := ■^^-{W''+^>(i/p).i°g</'(|i^.-iin+i°g(i/^^ obtaining J,,; we stop at iteration k if log(*^+i) il/p) < 

^ " logC-) (l/p)+log{l/e) & A.' F t> \ /i^j 

log (pdFk-il). However, this process involves several technical difficulties since (i) the values of p, e change 
over iterations, and (ii) the (weighted) measure becomes somewhat intricate, as it should consist of 3"^ and 
all sets of heavy objects collected over all iterations. The author has several initial bounds obtained for this 
process, and she plans to finalize these details in the full version of this paper. This will tighten the current 
(probably, suboptimal) bounds stated in Corollaries 2.8 and 2.9. 

Another interesting problem that we plan to study is whether the log (1/e) factor in the enumerator of 
our bound can be removed. Li et al. [17] obtained such an improvement, where they reduced the previously 
known bound 0(log {l/{pe)) /e'^p) in [14, 26] to 0(log {l/p)/e^p), which they showed to be optimal in 
the worst case. This improvement is derived by applying the chaining method of Kolmogorov. Roughly 
speaking, in this technique the standard union bound over a set of events (defined in some probability space) 
is replaced by a tighter bound when considering only a relatively small subset of events, each of which is 
"distinct" in some sense (such a subset is also called an "e-packing" [12]). It is a challenging open problem 
to combine our machinery with the chaining method. Specifically, does an e-packing exist in our scenario? 
If so, can one apply the Local Lemma of Lovasz on the corresponding events? 

Last but not least is the implications of our approach to the bounds on combinatorial discrepancy for 
well-behaved range spaces. In particular, even just the case of points and halfplane ranges is already chal- 
lenging. Har-Peled and Sharir [13] showed that in such range spaces (P, 'K) the discrepancy x{'t ^ P) of 
each range r G !K is only 0(|r n logn), where n = \P\. This property eventually yields the improved 
bound for relative (p, e) -approximations in this scenario; recall that the improvement in [13] is with respect 
to the dependency on e. If the factor log n in xi^ ^ P) can be reduced to o(log n) then this will yield an 
improvement in the parameter p as well. Nevertheless, we were unable to apply our technique on this setting 
so far, due to the differences in our probabilistic model and the one applied in [13]. 
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A The Construction 



A.l Proof of Lemma 2.2: 

For simplicity of presentation, we bound Proh{Ar}, for all r G So by ( iog(i/p)(/>( | 5|) ) ' which clearly 
holds by the bound given in the proof of Proposition 2. 1. Thus (3) holds for alH = 0, . . . , log 

Fix a range r and a layer Sj. We first observe that < :s.{Ar) < 1- The lower bound is trivial. For the 

upper bound, we obtain from (3) that x{Ar) < ( , — n/ u^flTh ) ' which is smaller than 1 since B > 1 



Jog{l/p)<t>{\J\)^ 

(and is chosen to be sufficiently large) and p < 1/8 by our assumption. 

We next consider all events Ar' with A^' ~ Ar, A^i / At-. Let i be the layer of the range r, and let 
i be the layer of r'. The corresponding product ^^j^^ ^ (1 — x(At-/)) in the right-hand side of (4) 
can be spelled out for r as 

log(l/p) 

n n (i-x(^r'))- 

j=0 A^,^AT,A^,^Ar 



A .en 



We lower bound separately the sub-products involving layers Sj with i < j and layers $j with i > j. 

(i) i < j: In this case, since r contains at most 0(Aj) light objects (as noted above, r may also contains 
at most 0(Aj) heavy objects, but they are ignored at that part of the analysis), each of which partici- 
pates in at most 0[4>{\'J\)/S.j'^'^'^) ranges r' of Sj, it follows that r has a degree-of-dependency at most 
a ■ (/)(|9'|)AjAj'^"''^, for some absolute constant a > 0, with the ranges of Sj. We thus obtain: 



n a-.,A.)),h- j^^^j^ 



A ,~Ar,A ,^A 



1 



We simplify the right-hand side using the inequality (1 — y)^ > 1 — zy, for < y < 1, and the fact that 
Afe = — ^°g(^/(P^)) ^ for each k = 0, . . . , log (l/p)- Specifically, we have 

a . <^(|:r|)A,A/+2 < a • <A(|:?I)A/+^ = a' • 0(1^1) ^g^^^^^""'^^ ' 
for another constant a'. We then have 

log(l/p),^(|3-|)J j - l^log(l/p)<A(|:J|)J 

^ ^ - l,iog(i/p)^(|:^|) J • « • ^(1^1) ^) • 

By "stealing" a small portion /3 • 2-'"^ from the exponent B ■ 2^^^, for some sufficiently large constant /3 > 0, 
we can assume that 

f ^'""e , ^,,^„ 2(-+3)nog^+3(i/(pg)) 
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By choosing B sufficiently large, we thus obtain: 

n (i-x(^.o)>i 

TeSi,T'eSj,j>i 



\og{l/p)cP{\J\) 



for a suitable constant < B' < B that can be made arbitrarily large by choosing B sufficiently large. 
Having such a choice for B', and recalling that p < 1/8, we obtain that the latter expression is greater than 

1 — (i)^ , where k > is a large constant, whose choice depends on B' (and thus on B). 

(ii) i > j: As in the previous case, r has a degree-of-dependency at most a ■ (/>(|3"|)Ai Aj'^"'"^ with the layers 
of Sj, and we have: 



n (i-x(A.o)> 1 



^2^-1 X a-0(|S-l)A,A/+2 



Choosing B sufficiently large, we can assume that ( iog(i/p)0^(|j-|) j ^ 1/2- We can now use the 

inequality (1 — t/Y > e"^^^, for < y < 1/2. Proceeding as in the previous case, we have 



a-0(|J|)A,Af2 = a'.<^(|J|) 



£2(c+3) ' 

1-1 . A I > 



so 



log(l/p)</.(|3-|) y 



Applying an exponent stealing similar to the one above, we can cancel most of the other part of the expres- 
sion, and end up with a lower bound of the form: 



n (1 - x(A,0) > exp \ -2' 



A^,~At,A^,^At 



iogii/p)<pim) 



where < B" < B is another constant sufficiently close to B. We can in fact assume that B" = B' by 
replacing the lai^ger of them by the smaller. Choosing B (and thus B') to be sufficiently large, and using, as 
above, the inequality p < 1/8, we obtain that 



n (l-x(^.0)>exp|-2'-Q] 



2-'- 



A^ir~.Ar,A^,^A 

where k > can be chosen to be the same constant as in case (i). 
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Putting the bounds together: We thus obtain, for a fixed range r, 

n n (i-x(A.))>n«p -2- (i) ■ n 

r'es,- 



For the first product we have 



i=o 

For the latter product we have 



Hexp I -2^ . I = exp | -2^^ (1^ ^ > exp{-2V(^- 1)} . 



log(l/p) / log(l/p) /]^\2^"^ 2^ 



n >- i E a 



This implies that 



> 1- ^.-^ > 1 



kJ - K^'-' - 1 - JE- 1 



n n (1 - ^(^-')) > 1 - 7^ • exp {-2V(v^ - 1)} . 

j=0 A^,^Ar,A^,j^A^ V V / 



Given that k is chosen to be sufficiently large we obtain 

log (1/p) 

exp {2'+^} ■ Proh{Ar} ■ JJ (1 - x(^^/)) > Prob{yi^}, (7) 

j=0 A^,r^AT,A^,T^Ar 

as asserted. This completes the proof of the lemma. 
A.2 Proof of Corollary 2.3 

To this end, we strengthen (5) to include the event 23 : |3"i| > (1 + 7) Exp{|Fi}, where 7 is a sufficiently 
large constant.^ That is, we show that there exists an additional assignment x' : S — (0, 1), such that 

(1 - Prob{S}) /\ i /\ (1 - Proh{Ar}) i > (1 - x'(B)) J] (1 - x(yiO) > 0, (8) 

implying that all complementary events occur with a positive probability. By construction, for each each 
event Ar € £, we have Ar ~ 25. We thus need to modify (4) as follows: 

Prob{^^} < x(yi^) • (1 - x'(!B)) ]J (1 - x(^^/)), for each Ar G £, and (9) 



*If fact, it is sufficient to show that Prob{B} is smaller than /\''!f,j*^^''^ Atgs ~ Prob{ylr}). However, in our analysis we 
include the event 3 into the local lemma in order to eventually be able to apply it in a constructive manner. 
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Prob{S} < x'(B) • n ~ x(^^)). 

Indeed, we obtain using Chemoff's bound once again, the fact that Exp{ 1 9"i I } = 7r|L|, and |L| > |3"|/2: 
Prob{S} = Prob{|Ji| > (l+7)Exp{|Fi|}} <exp{-7V3-Exp{|Ji|}} (10) 



< 



exp 



DY max{loglog log (0(| J|)} + log (1/e) 



6 e^p 
By (5) and the assignment of x(yiT-) in Lemma 2.2 we obtain: 

log (1/p) 

l\ Vvoh{Ar}) > n (1 - ^{-^-i) = n (1 - 2'^' • Prob{yi4) . 

i=0 T&i ylT6£ A-r&B. 

Using (3) and the fact that there are only 0(|3"|(/)(|3"|)A^) ranges in layer i, the latter term is lower bounded 



n 1 



i=0 



log(l/p)</<(|J|) j 



log (1/p) f / \ 



> TT ; o( e^'^'-e \ , 2- log'^+i 



j=0 



log(l/p)(/)(|J|) / e2(c+i)^ 



for two absolute constants a, a' > 0. Applying once again "exponent stealing" similarly to the proof of 
Lemma 2.2, we obtain that this term is at least 



exp 



^ [iog{i/p)cpm) ■pr-'^^pj I l^iog(i/p)0(|j|)J 



for a suitable constant < B' < B that can be made arbitrarily large by choosing B sufficiently large. 
Comparing the latter exponent with that of (10), we conclude that 

L>72 log log (1/p) + log (l/e) \ \ l^f e^/^ • ^ \ ^ 1 



exp <^ — • 5 } < exp < - 2 



e'p J 1 \" \log(l/p)<^(|J|) j p 



when 7 is chosen to be sufficiently large. We now put x'(!B) to be 

B'/2 



exp < 



2 



e^/^ • g \ l\ maxjloglog (1/p), log ((/'(|:J|)} + log (l/e) 

log(l/p)0(|J|) J ■p)"~6 



In other words, we set x'(23) to be the ratio between the upper bound on Prob{S} and the lower bound 
on OyireE ~ ^{-^r})- Due to this property we always have x'(15) • Hyi^ee (-*- ~ ^i-^r}) > ProblS}, 
and thus the second part of Inequality (9) is satisfied. We also note that when 7 is chosen to be sufficiently 
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large, we should have x'('B) < 1. In fact, x'(23) can be made arbitrarily small (by choosing 7 sufficiently 
large), and then we have, say, x'(23) < 1/2. This implies that we can modify Inequality (7) from the proof 
of Lemma 2.2 so that it now satisfies 

log (1/p) 

exp {2^+1} • Prob{A^} • (1 - x'(S)) J] 11 " ^(^^')) ^ Prob{A^}, 

which shows the first part of (9). 
A.3 Applications 

A.3.1 Points and Axis-Parallel Boxes in Two and Three Dimensions 

We begin with the two-dimensional case. It is well known that a set P of n points in the plane admits 
0(n^) rectangular "empty" ranges (that is, these ranges do not contain any point of P in their interior); see, 
e.g., [3]. Thus we resort to the technique presented in [3, 10] instead. We use the following property shown 
in [10] (and based on the analysis in [3]): 

Lemma A.l (Ene et al. [10]). Given a set P of n points in the plane and a parameter k > 0, one can 
compute a set Rk of 0{k'^n log n) axis-parallel rectangles (each of which is "anchored" either on its right 
side or its left side to a vertical line), such that for any axis-parallel rectangle r, if \r r\ P\ < k, then there 
exists two rectangles ri,r2 E Rk such that |(ri U H P| = |r n P|. 

Thus in this case c = 2 and (/){■) is the log(-) function, so the range space [P, Rn) is well-behaved. In 
order to bound the relative eiTor for the original (non-anchored) rectangular ranges we proceed as follows. 
Let r be such a range realized by a rectangle r, let ri, r2 be its two corresponding portions satisfying 
the property in Lemma A.l, and let ri, T2 be {ri n P}, {r2 n P}, respectively. We now replace P by 
the corresponding sample J, and, as before, denote by T the set of the resulting ranges, and then apply 
Theorem 2.4 in order to obtain: 

[l- e) - ep< r— < — — - — (1 + e) + ep, 



for i = 1, 2. By combining these inequalities for ti,T2, we obtain: 

IrnJL , IrnJil + 7r|Tni?| IrnJL 
' '-(1 -e) - 2ep < J ^ < ' ' (1 + g) + 2£p, 



|3-| ^ ' vrlJl - |J| 

This implies that 3"i U i7 is a relative (p, 3e) -approximation for (3", T) (in the above "weighted sense"), and 
1:^1 UH\= 0((loglog + log {l/e))/e^p). 

When P is a set of points in three dimensions, one can obtain similar properties to those in Lemma A.l, 
derived from the analysis in [3]. In this case one can compute a set Bk of 0{k'^n log^ n) axis-parallel boxes 
such that for any axis-parallel box 6, if |fe n P| < k, then there exists eight boxes 61, ... 63 G Bj^ such that 
|(U?^]^6j) n P| = |6 n P|. We omit the straightforward details in this version. 

Thus in this case c = 2 and is the log'^(-) function, so the range space (P, Bn) is well-behaved. As 
above, put ti := {hi n P}, i = 1, . . . , 8. Replacing P with 3" once again, applying Theorem 2.4, and then 
combining the resulting inequalities for Tj, i = 1, . . . , 8: 

' '-(1 - e) - Sep < ^J— p < ' ' (1 + e) + Sep, 



m ' ' vrlJl - |J| 
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and thus U is a relative (p, 9e) -approximation for (9", T), and its size has the same asymptotic bound 
as in the two-dimensional case. 



Points and a-fat triangles in the plane. When P is a set of n points in the plane and the ranges are a-fat 
triangles (that is, triangles for which each of their angles is at least a) the analysis in [10] implies: 

Lemma A.2 (Ene et al. [10]). Given a set Pofn points in the plane, a parameter k, and a constant a > 0, 
one can compute, in polynomial time, a set ofO{k^n log^ n) regions, such that for any a-fat triangle A, 
/fl A n P| < k, then there exists (at most) 9 regions in 7k whose union has the same intersection with P as 
A does. 

Using similar arguments and notation as in the case for axis-parallel boxes (here we have c = 3 and (/)(•) 
is the log^(-) function), we obtain: 

' -(1 - e) - 9ep < ^J— p < ' ' (1 + e) + 9ep, 



\J\ ' vrlJl - |J| 

and thus is 3"i U i7 is a relative (p, 10e)-approximation for (3", T), and its size has the same asymptotic 
bound as in the previous cases. 

We now appropriately scale the parameters e, p, for each of the above three settings, in order to obtain 
Corollary 2.7. 

A.3.2 Points and Halfspaces in Two and Three Dimensions 

Let P be a set of n points, and let Ji denote all halfspace ranges. Using the random sampling theory of 
Clai^kson and Shor [7], for any subset of m points of P, the number of halfspaces of size at most k, for any 
integer parameter < A; < m, is 0{mk + m) in two dimensions and 0{mk^ + m) in three dimensions. It 
thus follows that in both cases (P, "K) is a well-behaved range space, and we can thus apply Theorem 2.6 
for these cases, thereby showing Corollary 2.8. 

A.3.3 Planar Regions of Nearly-Linear Union Complexity and Points 

Let P be a set of n (closed) connected planar regions, and let U(P) = [J R denote the union of R. The 
combinatorial complexity of U(P) is the number of vertices and edges of the amngement A{R) that appeal" 
along dUlR). For r < n, let u(r) denote the maximum complexity of the union of any subset of r regions in 
R, measured as above. We assume that u(r) is nearly linear, i.e., u(r) < rtp{r), where tp{-) is a (sublinear) 
slowly growing function. 

We now consider the dual range space defined on R and all points in the plane. In fact, it is fairly 
standard to show that these ranges correspond to all faces in the arrangement A{R), where each range is 
the subset of regions covering a fixed face of A{R). Using a standard application of the Clarkson-Shor 
technique [7], for any subset of r regions of R, the number of such faces of size at most k, for any integer 
parameter > 0, is 0{k'^\i{r /k)), or 0{krip{r /k)). Applying Theorem 2.6, we obtain: 

Corollary A.3. Let R be a set of n planar regions such that the union complexity of any r of them is 
u(r) = r(p{r), where (/?(•) is a (sublinear) slowly growing function. Then any dual range space on R and a 
set of points in the plane admits a (weighted) relative (p, e)-approximation of size 

Q . ' max{ log log ( 1 /p) , log 99(log ( 1 / (ep) ) /g^p) } + log ( 1 /g) ) 

e^p 
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for any 0<p<l, 0<e<l. This relative approximation can be constructed in expected polynomial time. 

We now present several standard families with this property, state their union complexity, and then apply 
Corollary A.3 for each of these cases in order to conclude Corollary 2.9. 

Pseudo-disks. In a set of pseudo-disks, the boundary of any pair of regions are either disjoint or cross 
twice. In this case u(r) = 0{r) [18]. 

a-fat triangles. Recall that a triangle is a-fat if each of its angles is at least a. In this case u(n) = 
0(r log* r), where the constant of proportionality depends on a [2, 11]. When the triangles have roughly 
the same size, this complexity reduces to 0{r) [22]. When on the sides of these triangles is unbounded (in 
which case the triangles become a-fat wedges), this complexity also reduces to 0{r) [9]. 

Locally 7-fat objects. Given a parameter < 7 < 1, an object is locally j-fat if, for any disk D whose 
center lies in o, such that D does not fully contain o in its interior, we have area(D Ho) > 7 • area(L') , where 
D n o is the connected component of D n o that contains the center of D. We also assume that the boundary 
of each of the given objects has only 0(1) locally x-extreme points, and that the boundaries of any pair of 
objects intersect in at most s points, for some constant s. It is then shown in [2] that u(r) = 0{r ■ 2*^('°s* '')), 
where the constant of proportionality (in the linear term) depends on 7. 
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